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ULTRASONIC TRANSMITTING AND RECEIVING APPARATUS AND 
ULTRASONIC TRANSMITTING AND RECEIVING METHOD 

5 

BACKGROUND OF THE INVENTION 
Field of the invention 

The present invention relates to an ultrasonic 
transmitting and receiving apparatus and an ultrasonic 
10 transmitting and receiving method to be used for obtaining 
ultrasonic images by transmitting ultrasonic waves and 
receiving ultrasonic echoes . 

Description of a Related Art 

15 Conventionally , in order to obtain a two-dimensional 

or three-dimensional image, for example, like B -mode scanning, 
ultrasonic beams are transmitted one by one to scan an object 
to be inspected, and a two- dimensional or three-dimensional 
image is synthesized based on the obtained image information . 

20 However, according to such method, since time lag between 
frames is large, images in different times are synthesized 
and the synthesized i mage is blurred when imaging moving parts . 
Especially, when observing a part moving hard such as a 
circulatory organ, real time operation of at least 30 frames 

25 per second is required. In order to obtain ultrasonic images 
in real time, the imaging region within the object must be 
scanned at high speed, and the method of sequential scanning 
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with a single ultrasonic beam is too late for the requirement. 

As a countermeasure against such problem, a technology 
of transmitting and receiving plural ultrasonic beams 
simultaneously toward many directions from an ultrasonic 
5 transducer array in which plural ultrasonic transducers are 
arranged two-dimensionally is under study. 

Fig. 15A schematically shows a state in which ultrasonic 
beams are transmitted from an ultrasonic transducer array 
included in an ultrasonic probe that is generally used. An 

10 ultrasonic transducer array 100 is f abricatedby , for example , 
linearly arranging a number of ultrasonic transducers 101 . 
As an element (ultrasonic transducer) used for transmission 
and reception of ultrasonic waves, an element in which 
electrodes are formed on both ends of a material having a 

15 piezoelectric property (piezoelectric element) that includes 
piezoelectric ceramic represented by PZT (Pb (lead) zirconate 
titanate) or a macromolecule piezoelectric element 
represented by PVDF (polyvinylidene difluoride) has been 
generally used. 

20 Drive signal generating circuits including pulsers etc. 

are connected to these ultrasonic transducers 101, 
respectively. Applying a voltage to such ultrasonic 
transducer 101, the piezoelectric element expands and 
contracts by piezoelectric effect to generate ultrasonic 

25 waves . By driving plural ultrasonic transducers 101 at 

predetermined time intervals, spherical waves transmitted 
from the respective ultrasonic transducers 101 are 
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synthesized and the ultrasonic beam having a focal point formed 
in a desired direction and a desired depth can be transmitted. 

In addition, as shown in Fig. 15B, by applying two sets 
of timing pulses of A pulse and B pulse to one set of elements , 
5 an ultrasonic beam A and an ultrasonic beam B can be transmitted 
simultaneously in different directions. Note that, in the 
case where the A pulse and B pulse overlap, a common pulse 
as an addition result may be generated. 

Hereinafter, thus simultaneously transmitted and 
10 received plural ultrasonic beams are referred to as 
"multi-beam" . 

Now, in ultrasonic imaging, side lobes produced when 
transmitting ultrasonic beams become problematic. When an 
ultrasonic beam having directivity is transmitted, in spatial 

15 distribution of acoustic pressure intensity, a local maximum 
that occurs on the center axis in a transmitting direction 
is referred to as "main lobe" , and a local maximum that occurs 
in a direction other than that is referred to as "side lobe" . 
The side lobe is produced depending on the relationship between 

20 the element pitch of the ultrasonic transducers and ultrasonic 
frequency (such side lobe is referred to as "grating lobe") , 
or produced by the unwanted vibration of the ultrasonic 
transducers. Normally, an ultrasonic echo received by the 
ultrasonic transducer is subjected to signal processing as 

25 an echo that has propagated from the direction of the main 
lobe. On this account, in the case where the side lobe 
component is large or there is a strong reflector in the side 
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lobe direction, an artifact (virtual image) is generated, 
and the image quality of the ultrasonic image is degraded. 

In order to avoid such phenomenon, measures have been 
taken for suppressing the side lobes, such as improvement 
5 in the delay accuracy of the transmitted and reception beams 
andminiaturization of the element . However , there are limits 
to these techniques , and the side lobes have not been reduced 
to the sufficient level. Further, in the case where the 
ultrasonic beam is largely steered, or ultrasonic waves are 

10 subjected to multi-beam transmission, the level of the side 
lobe component becomes higher, and accordingly, it becomes 
more difficult to reduce it. Therefore, the effect on the 
image quality becomes a major problem. 

U.S. Patent No. 6,179,780 discloses the following 

15 technology for reducing the effect by the side lobes when 
performing multi-beam transmission and reception of 
ultrasonic waves. That is, a method of forming plural 
reception beams for one transmission beam, or a method of 
identifying transmission beams by changing frequencies of 

20 plural transmission beams or coding transmission beams by 
using Barker code, Golaycode, etc. to relate them with received 
ultrasonic echoes are cited. Further, since there is a region 
called as "null line" where the acoustic pressure becomes 
generally zero between a main lobe and a side lobe, a method 

25 of performing alignment of the main lobe of another ultrasonic 
beam in the region , a method of simply separating the intervals 
of transmission beams, and a method of shifting the center 
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frequencies of the 'transmission beams are also cited. 
Although these methods are for suppressing the occurrence 
of the side lobes , there Is a limit to the degree that the 
side lobes can be reduced. 
5 By the way, recently, as one of greatly developed 

technologies in the medical Imaging technology field, X-ray 
CT (computed tomography) can be cited. The X-ray CT is a 
technology of applying X-rays from plural directions to the 
object and generating the tomogram of the object based on 

10 the X-rays transmitted through the object. However, since 
X-ray signals transmitted through the object from different 
hundreds of directions are required in order to obtain an 
X-ray CT image with high image quality, a complex and accurate 
mechanism for rotating the X-ray source and X-ray detecting 

15 unit, a huge amount of time and a high performance data 
processing unit for processing an enormous amount of data, 
etc. are indispensable . Accordingly, the apparatus becomes 
large scaled and expensive, and that has discouraged the 
widespread in general use . Further , the exposed dose of the 

20 X-ray applied to an object to be inspected has been also 
problematic . 

In order to solve such problems, Japanese Patent 
Application Publication JP-A- 9-1 61041 discloses a CT 
apparatus capable of obtaining high quality images at high 
25 speed even if the number of projection directions is small. 
In this CT apparatus, a generalized inverse matrix is stored 
which has been calculated in advance by performing singular 
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value decomposition on the projection model matrix expressing 
in a matrix the relationship between predetermined plural 
different transmission paths and influence coefficients as 
the degree of effect on the projection value by respective 
5 pixels within a calculation space partitioned in a lattice 
form when a transmission ray passes through respective 
transmission paths, and then, the data for image display is 
generated by using the projection values on intersections 
of the X-rays projected from different directions (the 

10 partitioned calculation space in a lattice form) and the 
generalized inverse matrix. 

However, this technique is for estimating the signals 
on the lattice points, where these transmission paths 
intersect, from the projection values in the different plural 

15 transmission paths (integration values of the transmission 
paths) , and such technology can not be directly applied to 
the ultrasonic imaging for obtaining information on the object 
from the acoustic pressure intensity of the ultrasonic echoes 
reflected from the object. 

20 

SUMMARY OF THE INVENTION 

The present invention has been achieved in view of the 
above-described problems . An object of the present invention 
is to provide an ultrasonic transmitting and receiving 
25 apparatus and an ultrasonic transmitting and receiving method 
capable of obtaining high image quality ultrasonic images 
with reduced effect of side lobes . 
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In order -bo solve the above-described problems, an 
ultrasonic transmitting and receiving apparatus according 
to the present Invention comprises : an ultrasonic transducer 
array Including plural ultrasonic transducers for 
5 transmitting ultrasonic waves and receiving ultrasonic echoes 
reflected from an object to be Inspected; drive signal 
generating means for generating drive signals for 
respectively driving the plural ultrasonic transducers; 
transmission control means for controlling the drive signal 

10 generating means such that ultrasonic waves transmitted from 
the plural ultrasonic transducers form at least one ultrasonic 
beam; signal processing means for performing reception 
focusing processing on plural detection signals obtained 
based on the ultrasonic echoes received by the plural 

15 ultrasonic transducers so as to form a reception focal point 
in at least one region within the object thereby obtaining 
plural detection signals relating to the at least one region; 
storage means for storing plural different acoustic pressure 
intensity prof lies set based on acoustic pressure intensity 

20 distribution formed by transmitting ultrasonic beams and the 
reception focusing processing to be performed in the signal 
processing means ; and calculating means for calculating image 
data relating to the at least one region on the basis of the 
plural detection signals relating to the at least one region 

25 and the plural different acoustic pressure intensity 
profiles . 

Here, the signal processing means may obtain plural 
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detection signals respectively relating to plural regions 
within the object, and the calculating means may calculate 
image data relating to the plural regions by using the plural 
detection signals and plural different acoustic pressure 
5 intensity profiles relating to the plural regions . Further, 
the signal processing means may obtain plural detection 
signals relating to one region within the object, and the 
calculating means may calculate image data relating to the 
one region by using the plural detection signals and plural 

10 different acoustic pressure intensity profiles relating to 
the one region . Furthermore , the signal processing means may 
obtain plural detection signals relating to each of plural 
regions within the object, and the calculating means may 
calculate image data relating to the plural regions by using 

15 the plural detection signals and the plural different acoustic 
pressure intensity profiles relating to the plural regions . 

An ultrasonic transmitting and receiving method 
according to the present invention is an ultrasonic 
transmitting and receiving method of transmitting and 

20 receiving ultrasonic waves using an ultrasonic transducer 
array including plural ultrasonic transducers for 
transmitting ultrasonic waves and receiving ultrasonic echoes 
reflected from an object to be inspected, and the method 
comprises the steps of : (a) transmitting at least one 

25 ultrasonic beam by driving the plural ultrasonic transducers , 
respectively; (b) obtaining plural detection signals relating 
to the at least one region within the object by performing 



reception focusing processing on the plural detection signals 
obtained based on the ultrasonic echoes received by the plural 
ultrasonic transducers so as to form a reception focal point 
in the at least one region ; (c) calculating image data relating 
5 to at least one region based on plural different acoustic 
pressure intensity profiles set based on acoustic pressure 
intensity distribution formed by transmitting ultrasonic 
beams and the reception focusing processing performed at step 
(b) and the plural detection signals relating to the at least 

10 one region. 

According to the present invention, since ultrasonic 
echo signals with no side lobe are estimated based on the 
acoustic pressure intensity profiles and the detection 
signals, high quality ultrasonic images with reduced effect 

15 of side lobes can be obtained. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram showing a constitution of an 
ultrasonic transmitting and receiving apparatus according 
20 to the first embodiment of the present invention ; 

Fig. 2A is a schematic diagram showing a state in which 
plural ultrasonic beams are transmitted from an ultrasonic 
transducer, and Fig. 2B shows the profiles of transmission 
beams formed within the object by transmitting plural 
25 ultrasonic beams; 

Fig. 3A is a schematic diagram showing plural reception 
beams having reception focal points formed therein, and Fig 
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3B shows profiles of -the reception beams; 

Fig. 4 shows acoustic pressure Intensity profiles 
obtained based on the profiles of the transmission beams and 
the profiles of the reception beams; 
5 Figs . 5Aand5B show acoustic pressure Intensity prof lies 

obtained by transmitting ultrasonic beams In three directions 
and performing reception focusing processing with respect 
to five directions; 

Fig. 6 show acoustic pressure intensity profiles 
10 obtained by performing reception focusing processing with 
respect to each of three regions; 

Fig. 7 is a flowchart showing an ultrasonic transmitting 
and receiving method according to the first embodiment of 
the present invention; 
15 Fig. 8 is a block diagram showing an ultrasonic 

transmitting and receiving apparatus according to the second 
embodiment of the present invention; 

Fig. 9 is a flowchart showing an ultrasonic transmitting 
and receiving method according to the second embodiment of 
20 the present invention; 

Figs . 10A and 10B are diagrams for explanation of an 
ultrasonic transmitting and receiving method according to 
the third embodiment of the present invention ; 

Figs . 11A and 11B are diagrams for explanation of an 
25 ultrasonic transmitting and receiving method according to 
the fourth embodiment of the present invention; 

Fig. 12 is a block diagram showing a modified example 
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of the ultrasonic transmitting and receiving apparatus shown 
in Fig. 8; 

Fig . 13 is a diagram for explanation of a modified example 
of the ultrasonic transmitting and receiving methods in the 
5 first to fifth embodiments of the present invention; 

Fig. 14 is a flowchart showing the modified example of 
the ultrasonic transmitting and receiving methods in the first 
to fifth embodiments of the present invention; and 

Fig. 15A is a diagram showing a state in which ultrasonic 
10 beams are transmitted from an ultrasonic transducer array 
included in a conventional ultrasonic probe, and Fig. 15B 
shows timing pulses applied to the respective ultrasonic 
transducers . 

15 DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Hereinafter, embodiments of the present invention will 
be described in detail by referring to the drawings . The same 
component elements are assigned with the same reference 
numerals and the descriptions thereof will be omitted. 

20 Fig. 1 is a block diagram showing a constitution of an 

ultrasonic transmitting and receiving apparatus according 
to the first embodiment of the present invention. 

The ultrasonic transmitting and receiving apparatus 
according to this embodiment includes an ultrasonic 

25 transducer array 10. The ultrasonic transducer array 10 
includes plural ultrasonic transducers arranged in a matrix 
form, for example. The ultrasonic transducer array 10 
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transmits ultrasonic beams based on driving signals applied 
thereto , and receives propagating ultrasonic waves to output 
detection signals . Each ultrasonic transducer is 
constituted by a vibrator in which electrodes are formed on 
5 both ends of a material having a piezoelectric property 
(piezoelectric element) that includes piezoelectric ceramic 
represented by PZT (Pb (lead) zirconate titanate) or a 
macromolecule piezoelectric element represented by PVDF 
(polyvinylidene di fluoride) , for example. Applying a 

10 voltage of a pulsed electrical signal or continuous wave 
electrical signal to the electrodes of such vibrator, the 
piezoelectric element expands and contracts. By the 
expansion and contraction, pulsed or continuous ultrasonic 
waves are generated from respective vibrators, and these 

15 ultrasonic waves are synthesized to form an ultrasonic beam. 
Further, the vibrators expand and contract by receiving 
propagating ultrasonic waves and generate electrical signals . 
These electrical signals are outputted as detection signals 
of ultrasonic waves . 

20 Alternatively, as the ultrasonic transducers, plural 

kinds of elements according to different ultrasonic 
conversion methods may be used. For example, the 
above-described vibrator is used as an element for 
transmitting ultrasonic waves, and a photo-detection type 

25 ultrasonic transducer is used as an element for receiving 
ultrasonic waves . The photo-detection type ultrasonic 
transducer is for detecting an ultrasonic signal by converting 
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it into an optical signal, and, for example, constituted by 
a Fabry-Perot resonator or fiber Bragg grating. 

Further, the ultrasonic transmitting and receiving 
apparatus includes a scanning control unit 11, a transmission 
5 delay pattern storage unit 12 , a transmission control unit 
13, a drive signal generating unit 14, and a transmission 
and reception switching unit 15 . The scanning control unit 
11 sequentially sets transmitting directions of ultrasonic 
beams and receiving directions of ultrasonic echoes . The 

10 transmission delay pattern storage unit 12 has stored plural 
transmission delay patterns to be used when the ultrasonic 
beams are formed. The transmission control unit 13 selects 
a transmission delay pattern from among plural transmission 
delay patterns stored in the transmission delay pattern 

15 storage unit 12 according to the transmitting directions set 
in the scanning control unit 11, and sets the delay times 
to be respectively provided to the plural ultrasonic 
transducers on the basis of the selected transmission delay 
pattern . 

20 The drive signal generating unit 14 is constituted by, 

for example, plural pulsers corresponding to the plural 
ultrasonic transducers . Each of the plural pulsers generates 
a drive signal based on the delay time set in the transmission 
control unit 13 . The transmission and reception switching 

25 unit 15 switches the generation of drive signals in the drive 
signal generating unit 14 and capture of the detection signals 
in a signal processing unit 21 at timing according to the 
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control of -the scanning control unit 11 . By restricting the 
time periods for reading detection signals as described above , 
echo signals reflected from a particular depth of the object 
can be detected. 
5 Further, the ultrasonic transmitting and receiving 

apparatus includes the signal processing unit 21, a primary 
storage unit 22, a reception delay pattern storage unit 23, 
a reception control unit 24, a secondary storage unit 25, 
an intensity profile storage unit 31, an essential data 

10 calculating unit 32, a tertiary storage unit 33, an image 
processing unit 34, and a display unit 35. 

The signal processing unit 21 includes plural channels 
corresponding to the plural ultrasonic transducers . Each of 
the plural channels of the signal processing unit 21 captures 

15 the detection signal outputted from the corresponding 

ultrasonic transducer with predetermined timing, andperf onus 
signal processing such as logarithmic amplification, 
demodulation, STC (sensitivity time control) , filtering 
processing, A/D conversion. The primary storage unit 22 

20 stores detection signals subjected to signal processing in 
the signal processing unit 21 with respect to each channel 
in chronological order. 

The reception delay pattern storage unit 23 has stored 
plural reception delay patterns to be used when reception 

25 focusing processing is performed on the plural detection 
signals outputted from the plural ultrasonic transducers . 
The reception control unit 24 selects a reception delay pattern 
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from plural reception delay patterns stored in the reception 
delay pattern storage unit 23 according to the reception 
direction set in the scanning control unit 11 , and performs 
reception focusing processing by providing delays to the 
5 plural detection signals based on the reception delay pattern 
and adding them to each other. By the reception focusing 
processing, sound ray data in which the focal point of the 
ultrasonic echo is narrowed down is formed. The secondary 
storage unit 25 stores the sound ray data formed in the 

10 reception control unit 24. 

The intensity profile storage unit 31 has stored plural 
different acoustic pressure intensity profiles, which will 
be also simply referred to as "intensity profiles" or 
"profiles" . The intensity profiles represent acoustic 

15 pressure intensity or acoustic pressure intensity ratios in 
plural regions included in the surface where the ultrasonic 
beams reach in a predetermined time elapsed from being 
tr an smi t ted . 

Hereinafter, the surface where the ultrasonic beams 
20 reach in predetermined time elapsed from being transmitted 
is referred to as "isochronal surface" . It is conceivable 
that, in plural regions included in one isochronal surface, 
distances from the transmitting position of the ultrasonic 
beam are substantially equal. The intensity profiles are 
25 expressed by the function of the positions (directions) of 
the plural regions included in the isochronal surface and 
the acoustic pressure intensity or acoustic pressure 
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intensity ratios therein. 

It is assumed that, in the space in which ultrasonic 
waves are transmitted and received, spatial distribution of 
acoustic pressure intensity (acoustic field) is formed by 
the acoustic pressure intensity distribution formed by 
transmitting ultrasonic beams and the reception focusing 
processing performed in the reception control unit 24. In 
this embodiment, such spatial distribution of acoustic 
pressure intensity is utilized in order to eliminate the effect 
of side lobes in the detection signals . 

The plural intensity profiles stored in the intensity 
profile storage unit 31 are set based on (i) the transmitting 
directions of ultrasonic beams and (ii) the reception focusing 
processing of ultrasonic echoes, that is, receiving 
directions and depths of the reception focal points . Further , 
as factors for setting the acoustic pressure intensity prof ile , 
a number of used elements from among the elements included 
in the ultrasonic transducer array 10, an element pitch, an 
aperture diameter of the used elements , or aperture conditions 
including a weighting factor within the aperture may be 
included. Such intensity profile can be obtained by 
simulating an acoustic field based on the above-described 
aperture conditions, transmission conditions including the 
transmission delay pattern, and reception conditions 
including the reception delay pattern. On the other hand, 
the intensity profile may be obtained by transmitting and 
receiving ultrasonic beams to the scattering phantom based 
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on these conditions and using intensity ratio of ultrasonic 
echoes obtained by the measurement. 

The essential data calculating unit 32 selects a profile 
from the plural profiles stored in the intensity profile 
5 storage unit 31 based on the transmitting direction and 
receiving direction set in the scanning control unit 11 , and 
calculates image data with the suppressed side lobe components 
(essential data) based on the profile and measurement data 
corresponding to a selected isochronal surface from among 

10 the sound data stored in the secondary storage unit 25. The 
tertiary storage unit 33 stores the image data calculated 
in the essential data calculating unit 32 . 

The image processing unit 34 constructs two-dimensional 
or three-dimensional image data based on the image data stored 

15 in the tertiary storage unit 33 , and performs image processing 
such as gain adjustment, contrast adjustment, gradation 
processing, response enhancement processing, interpolation 
processing, etc. and further performs scanning conversion 
on the image data. The display unit 35 displays ultrasonic 

20 images based on the image data processed in the image processing 
unit 34. The display unit 35 includes a display device such 
as a CRT or an LCD, for example. 

Next, the principle of the ultrasonic transmission and 
reception in this embodiment will be described by referring 

25 to Figs. 2A to 4. 

As shown in Fig. 2A, acquisition of essential data x A/ 
x B , x c relating to regions A, B, C on the isochronal surface 
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PHI in TA direction, TB direction, TC direction will be 
discussed. Here, the isochronal surfaces PHI, PH2 , PH3 are 
surfaces where the ultrasonic beams reach in time periods 
tl , t2 , t3 f rombeing transmitted, respectively . In addition , 
5 the essential data is ideal data representing information 
on the respective regions, in which the effect by the side 
lobe component is reduced. When ultrasonic beams are 
subjected to multi-beam transmission towards these three 
directions, at least three lobes (local maximums) occur in 

10 the acoustic field. Accordingly, the transmission beams are 
formed and transmitted so that strong lobes may be produced 
in TA direction, TB direction, TC direction shown in Fig. 
2A. In order to form desired transmission beams, a number 
of elements to be used from among the elements included in 

15 the ultrasonic transducer array 10, the element pitch, the 
aperture condition including an aperture diameter of the 
element to be used and weighting within the aperture, and 
the transmission condition including the transmission delay 
pattern may be used. Thereby, the acoustic field having the 

20 transmission beam profile as shown in Fig. 2B is formed in 
the space within the object. 

On the other hand, as shown in Fig. 3A, when ultrasonic 
echoes are received, reception focusing processing is 
performed on the detection signals so that reception focal 

25 points may be formed in a predetermined depth (for example, 
isochronal surface PHI) in the respective directions. Fig. 
3B shows profiles of reception beams RX1, RX2 , RX3 obtained 



by performing reception focusing processing with respect to 
each of TA direction, TB direction, TC direction. 

Accordingly, as shown in Fig. 4, intensity profiles 
obtained by actually transmitting ultrasonic beams and 
5 performing reception focusing on the received echoes are 
multiplication of the transmission beam profiles and 
reception beam profiles. In Fig. 4, ct A , Pa, Ya represent 
acoustic pressure intensity ratios corresponding to three 
lobes determined by the transmission beam TX and the reception 
10 beam RX1 . Similarly, a B , P B/ YBorac, p c , y c represent acoustic 
pressure intensity ratios corresponding three lobes 
determined by the transmission beam TX and the reception beam 
RX2 or RX3 . 

As shown in Fig. 4, even when the reception focusing 
15 processing is performed toward the predetermined direction, 
signal components from different directions are mixed in the 
measurement data obtained based on the actually received 
ultrasonic echoes. For example, the measurement data y A 
obtained based on the ultrasonic echoes received after time 
20 period ti x 2 from the transmission time includes a part of 
components of signals x B and x c respectively received from 
TB direction and TC direction other than the essential data 
x A representing information on the region A. Accordingly, 
assuming that the measurement data y A is obtained by performing 
25 weighted addition on the ultrasonic echo signals from these 
three directions (TA, TB, TC) with the above -de scribed 
acoustic pressure intensity ratios, the measurement data y A 
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Is expressed by the following equation. 

Ya = a A -x A + Pa'Xq + Ya-Xc 
Similarly, expressing measurement data y B/ yc by using the 
intensity profiles , the following simultaneous linear 
equations can be obtained. 

Ya = a A x A + p A x B + y A xc 

Yb = a B x A + Pb-Xb + y B -x c 

y c = ac*x A + p c x B + y c x c ... (1) 

Accordingly, by solving the simultaneous equations (1) 
by using the acoustic pressure intensity ratios a A , a B/ ... 
in the intensity profiles obtained in advance and the 
measurement data y A , y B/ y c , the essential data x A , x B , Xc can 
be obtained . 

Further, by performing such calculation with respect 
to all of the isochronal surfaces PHI , PH2 , . . . while varying 
transmitting directions TA, TB, ... of ultrasonic beams, 
essential data can be estimated with respect to all of the 
regions within the space. 

As above, the estimation of three essential data is 
explained by assuming that the ultrasonic beams are 
transmitted in three directions and the ultrasonic echo 
components propagating from different directions are mixed 
with three pieces of measurement data. However, the number 
of transmission beams, the number of measurement data to be 
used and the number of estimable essential data are not 
necessarily required to be the same. 

For example, as shown in Fig. 5A, five pieces of 
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measurement data may be acquired by performing multi-beam 
transmission toward the three directions of TA, TB, TC and 
performing reception focusing processing with respect to five 
directions of TA, TB, TC, TD, TE. In this case, as shown in 
5 Fig. 5B, five intensity prof iles are used. Accordingly, the 
five pieces of measurement data y A , y B , Yc, Yd, Yb in the 
isochronal surface PHI are expressed by the following 
equations by using the intensity profiles and essential data 
x A , x B , x c , x D , x E in the isochronal surface PHI. 
10 y A = ct A -x A + p A -x B + y A x c + 5 A -x D + e A -x B 

y B = <x B x A + p B x B + y B -x c + Se-x d + e B -x E 

yc = otc-x A + p c x B + y c -x c + 6 c -x D + e c -x E 

y D = ccd-x a + p D -x B + Yd-x c + 5d-x d + e D -x E 

y E = a E -x A + p E -x B + y E Xc + 8e-x d + e E x E ... (2) 
15 By solving the simultaneous equations (2) , the essential data 
x A to x E can be estimated. 

As described above, the estimation accuracy of the 
essential data can be raised by assuming that larger number 
of ultrasonic echo components are mixed in the measurement 
20 data and increasing the number of equations. For example, 
even in the case of performing multi-beam transmission in 
three directions, simultaneous equations may be constructed 
by assuming that ultrasonic echo components propagating from 
ten or more directions are mixed therein. 
25 By the way, in the simultaneous equations (1) and (2) , 

since the number of unknowns (essential data) and the number 
of equations are the same , the unknowns are uniquely determined . 
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However, in the case where at least one data having lower 
reliability is included in the plural pieces of measurement 
data, the estimation accuracy of all pieces of the essential 
data becomes lower. Accordingly, by constructing 
5 simultaneous equations so that the number of equations becomes 
larger than the number of unknowns , the stability of solution 
can be raised. 

In order to increase the number of equations constituting 
simultaneous equations, it is necessary to prepare plural 

10 pieces of measurement data relating to a predetermined region 
and plural different intensity profiles corresponding to 
their acquisition conditions. For that purpose, ultrasonic 
beams may be transmitted and received by varying one or all 
of the aperture condition of the ultrasonic transducer array 

15 10, transmission condition or reception condition of the 
ultrasonic beams. For example, by performing reception 
focusing processing on a series of detection signals with 
respect to the region A in the isochronal surface PHI in TA 
direction, which is shown in Fig. 2, by using different two 

20 kinds of reception delay patterns, two kinds of measurement 
data Yai and can be obtained. This is possible by setting 
reception focusing while varying the focal depth , for example . 
In addition, as shown in Fig. 6, two intensity profiles PR1 
(acoustic pressure intensity ratios a^, Pai/ Yai) and PR2 

25 (acousticpressure intensity ratios aA2/ 0a2/ Ya2> corresponding 
to the two kinds of reception delay patterns can be obtained. 
Similarly, plural pieces of measurement data and plural 
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intensity profiles are also obtained with respect to each 
of TB direction and TC direction . By expressing these pieces 
of measurement data by using essential data x A , x B , x c relating 
to the regions A, B, C and the intensity profiles , simultaneous 
5 equations (3) constituted by six equations can be constructed 
with respect to three unknowns . 
Yai = ctAi-x A + Pai-Xb + Yai-x c 

YA2 = CtA2-X A + Pa2-X B + Ya2'X c 

Ybi = a B i-x A + PbiXb + Ybi-Xc 
10 y B2 = a B2 x A + p B 2-x B + Yb2-x c 

Yci = Oci-Xa + PciXb + YciXc 
YC2 = Ctc2*X A + Pc2-x B + Yc2-x c . - . (3) 

Next, an ultrasonic transmitting and receiving method 

according to the first embodiment of the present invention 
15 will be described by referring to Figs . 1-4 and Fig. 7 . Fig. 

7 is a flowchart showing the transmitting and receiving method 

according to this embodiment. 

First, at step SI, the scanning control unit 11 sets 

transmitting directions of ultrasonic beams. Accordingly, 
20 the transmission control unit 13 sets delay times in the drive 

signal generating unit 14 based on the transmission delay 

patterns in response to the set transmitting directions as 

shown in Fig. 2B, for example. 

At step S2 , when the drive signal generating unit 14 
25 generates drive signals, the set number of ultrasonic beams 

are transmitted toward the set directions as shown in Fig. 

2A, for example. 
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At step S3, the switch is turned in the transmission 
and reception switching unit 15 , and the ultrasonic transducer 
array 10 receives ultrasonic echoes with predetermined timing 
Each of the plural ultrasonic transducers included in the 
5 ultrasonic transducer array 10 generates an electrical signal 
(detection signal) based on the received ultrasonic echoes. 

At step S4 , the signal processing unit 21 performs signal 
processing such as logarithmic amplification, STC, filtering 
processing, A/D conversion on the detection signal outputted 

10 from each of the plural ultrasonic transducers. At step S5, 
the signal processed detection signals (digital data) are 
sequentially stored in the primary storage unit 22 with respect 
to each line. 

At step S6, the reception control unit 24 performs 

15 reception focusing processing by using detection signals 
s tor ed in the primary s tor age uni 1 2 2 based on the predetermined 
reception delay patterns as shown in Fig. 3B, for example. 
Thereby, sound ray data corresponding to the transmitting 
directions set by the scanning control unit 11 are formed. 

20 The sound ray data formed at step S6 is stored in the secondary 
storage unit 25 (step S7) . 

At step S8 , the essential data calculating unit 32 
calculates essential data in a predetermined isochronal 
surface by using measurement data in the predetermined 

25 isochronal surface included in the sound ray data stored in 
the secondary storage unit 25, and intensity profiles (for 
example, see Fig. 4) stored in the intensity prof ile storage 
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unit 31 . 

The calculation of the essential data is performed as 
follows . The simultaneous equations described in the 
equations (1) to (3) can be expressed as the matrix relational 
expression (linear equation) as shown by the equation (4) . 

b = Ax ... (4) 

Where, "b" is a measurement data vector having measurement 
data y A , Yb, ... as components, and "x" is an essential data 
vector having essential data x A , x B , ... as components. 
Further, "A" is a matrix having "m" rows and "n" columns (m 
> n) and acoustic pressure intensity ratios a A , a B , ... in 
the intensity profiles as components. Specifically, the 
simultaneous equations (3) are expressed as follows: 
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In the case where m = n in the equation (4) , that is, 
the matrix "A" is a square matrix and a regular matrix, the 
vector "x" can be obtained by obtaining an inverse matrix 
A" 1 , and applying it from the left side of both sides of the 
equation (4) . In this case, as a method of obtaining the 
inverse matrix A" 1 , for example, the exact method including 
the sweeping-out method (Gaussian elimination method) and the 
LU decomposition method using it, or the iterative method 
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including the SOR method, etc. can be used. 

In addition, in the case where m > n in the equation 
(4) , by using the least-squares method, a plausible solution 
of the vector "x" can be obtained. 
5 The essential data calculated at step S8 is stored in 

the tertiary storage unit 33 (step S9) . 

At step S10, the image processing unit 34 constructs 
two-dimensional or three-dimensional image data based on the 
calculated essential data and performs image processing such 
10 as gain adjustment and gradation processing and scanning 
conversion. Then, at step Sll, the display unit 35 displays 
ultrasonic images on a display based on the image processed 
image data. 

As described above, in this embodiment, the side lobe 
15 component is eliminated from measurement data not by 
suppressing the occurrence of the side lobes in the 
transmission beams , but by estimating the essential data based 
on the measurement data and the acoustic pressure intensity 
profiles. Therefore, in the case where ultrasonic waves are 
20 subjected to multi-beam transmission, ultrasonic images with 
high image quality and reduced effect of side lobes can be 
obtained . 

Next , an ultrasonic transmitting and receiving apparatus 
according to the second embodiment of the present invention 
25 will be described by referring to Fig. 8. The ultrasonic 
transmitting and receiving apparatus according to this 
embodiment has a scanning part information input unit 40 and 
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a threshold setting unit 41. Further, the apparatus has an 
essential data calculating unit 42 instead of the essential 
data calculating unit 32 in Fig. 1. The rest of the 
constitution is the same as that shown in Fig. 1. 
5 The ultrasonic transmitting and receiving apparatus 

according to this embodiment uses the Moore-Penrose (MP) 
generalized inverse matrix when solving the simultaneous 
equations described in the equations (4) . Simultaneously, 
a matrix, a rank of which is reduced by discarding singular 

10 values, may be used. The MP generalized inverse matrix will 
be described later in detail . 

The scanning part information input unit 40 is provided 
on the operation panel, for example, and used when inputting 
a part of the object of ultrasonic imaging. The threshold 

15 setting unit 41 sets a threshold to be used when discarding 
singular values based on the scanning part information 
inputted from the scanning part information input unit 40. 
The threshold setting unit 41 may include a storage unit for 
storing plural thresholds corresponding plural parts of the 

20 object. 

Next, an ultrasonic transmitting and receiving method 
according to the second embodiment of the present invention 
will be described. Fig. 9 is a flowchart showing the 
ultrasonic transmitting and receiving method according to 
25 this embodiment. 

First, at step S21, the operator inputs scanning part 
information representing a part of the object of imaging by 
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using the scanning part information input unit 4 0 . In response 
to this, the threshold setting unit 41 sets a threshold to 
be used when discarding singular values based on the scanning 
part information inputted from the scanning part information 
5 input unit 40. 

Then, at steps SI to S7 , transmission and reception of 
ultrasonic waves are performed and the measurement data on 
a predetermined isochronal surface is stored in the secondary 
storage unit 25. The operation at steps SI to S7 is the same 
10 as described in the first embodiment of the present invention . 

Then, at step S22 , the essential data calculating unit 
42 calculates an essential data vector by using a matrix having 
acoustic pressure intensity ratios a A , a B , ... in intensity 
profiles as components, and a measurement data vector "b" 
15 stored in the secondary storage unit. 

Here, the method of calculating essential data at step 
S22 will be described in detail . 

In the special case where the matrix "A" is square and 
regular in the matrix relational expression (linear equation) 
20 shown by the equations (4) , as described in the first embodiment , 
the equation can be solved by directly obtaining the inverse 
matrix A" 1 . 

However, generally, the inverse matrix of the matrix 
"A" is not necessarily obtained uniquely . Especially, in the 
25 case where the number of equations is made larger than the 
number of unknowns (essential data) in order to raise the 
reliability of the solution as in the equations (3' ) , since 



the equations has extra conditions, the matrix "A" having 
acoustic pressure intensity ratios as components becomes a 
non-square matrix, and no inverse matrix exists. 

Here, a matrix "G" that satisfies AGA = A for the matrix 
5 "A" of m x n form is referred to as "generalized inverse matrix" . 
The generalized inverse matrix "G" gives least-squares 
solution to the essential data vector "x" of the equation 

(4) . However, the solution given by the generalized inverse 
matrix "G" is a general solution , and not determined uniquely . 

10 Accordingly, the following three conditions are given 

to the generalized inverse matrix "G" in addition to the 
condition of the generalized inverse matrix AGA = A. 
GAG = G (5) 
(GA) T = GA ... (6) 

15 (AG) T = AG ... (7) 

Where the matrix X T represents a transposedmatrix of the matrix 
"X" . 

The matrix "G" that satisfies all of these conditions 

(5) to (7) are referred to as "Moore-Penrose generalized 
20 inverse matrix", and represented by A + . The MP generalized 

inverse matrix A + is for uniquely determining generalized 
inverse matrix by using the condition that the norm of the 
vector "x" is minimum in the case where the least- square 
solution of the matrix "A" is indefinite. That is, by using 
25 the MP generalized inverse matrix, a plausible solution of 
simultaneous equations can be obtained. 

The MP generalized inverse matrix A + can be obtained by 
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the singular value decomposition (SVD) method. First, the 
m x n form matrix "A" is doubly diagonal! zed by the Householder 
transformation, and further diagonalized by the QR method 
to obtain r eigen values of X lr X 2 , - . . , A, r . Thereby, the 
5 singular value a± = is obtained. The r-th order diagonal 
matrix having this <J± as diagonal element is assumed as Z. 
Further, the m x r matrix or the n x r matrix formed by 
sequentially aligning normalized eigen vectors of 
corresponding symmetric matrix AA T and A T A are represented 
10 as "U", "V", respectively. By using these, the matrix "A" 
is decomposed as the equation (8) . 

A = UEV T ... (8) 
Here, the MP generalized inverse matrix A + is expressed by 
the following equation (9) . 
15 A + = VZ + U T ... (9) 

Where E + is a n x m diagonal matrix, and, assuming that the 
diagonal element of X is a±, the diagonal element a± + of E + 
is represented as follows . 
CTi + = l/<Ji (ai * 0) 
20 a ± + = 0 (di = 0) 

From the equations (4) and (9) , the essential data vector 
x is obtained by the equation (10) . 
x = A*b = VE + U T b . . . (10) 

In the actual calculation, it is desired that the rank 
25 of the matrix A is reduced, and the MP generalized inverse 
matrix is obtained by using the compressed matrix Aq, for 
relaxation of unstable conditions in simultaneous equations , 
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-that: is, for suppressing the noise Included In terms. 

Next, relaxation of unstable conditions In simultaneous 
equations will be described In detail . 

It is assumed that a vector with true measurement data 
5 including no error as components is b 0 , a vector having true 
essential data as components is x 0 , and a real blur function 
including no noise is Ao . The true measurement data vector 
b 0 is obtained by applying blur function Ao to the true essential 
data vector x 0 - 
10 b 0 = Aox 0 - . - (4' ) 

Correspondingly, the measurement data vector "b" and 
the blur function "A" in the equation (4) are expressed as 
follows . 

b = b 0 + k ... (11) 
15 A = Ao + 8A . . . (12) 

Where "k" is a vector having measurement data error as 
components. Further, SA is error of a blur function, and Ao 
is a real blur function . 

Using equations (4) , (4' ) , (11) , (12) , the measurement 
20 data "b" is expressed by the equation (13) . 
b = b 0 + k = Aox 0 + k 
= (A - 8A)x 0 + k 
= Ax 0 - SAxo + k 
b = Ax 0 + e . . . (13) 
25 where e is equivalent noise, and expressed as e = -SAx 0 + k. 

As described above, the linear inverse problem for 
obtaining an original image from a blurred image is substituted 
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with the estimation of the real essential data vector x 0 from 
the measurement data vector "b" . 

Here, the MP generalized inverse matrix A + is applied 
to both sides of the equation (13) . 
5 x = A*b = A + Ax 0 + A + e . . . (14) 

As clearly seen from the equation (14) , for example, in the 
case where, if "A" is slightly varied, "x" changes largely, 
that is, in the case where "A" is an adverse condition, A + e 
becomes extremely larger. In such case, it is not suitable 
10 to use x = A*b = x 0 + A + e as an estimated value. Accordingly, 
in order to obtain a solution with high accuracy, it is desired 
to reduce the rank of the matrix. It will be described later 
in detail that A is an adverse condition. 

The singular value decomposition of the matrices "A" 
15 and A + described in the equations (8) and (9) is represented 
as follows . 

A = ^ Kr w r vJ ...(15) 

r=l 

A + =i»X -(16) 
Where R represents a rank of the matrix "A" (R = Rank (A) ) , 

20 Ki, k 2 , . . . , k r , are singular values of the matrix "A" , and 
Ki > k 2 > ... k r > 0 . In addition, w r is a normalized eigen 
vector of corresponding AA T , v r is a normalized eigen vector 
of corresponding A T A. Further, K r _1 is an inverse number of 
the singular value . 

25 In the equation (15) , the matrix obtained by leaving 

larger Q values from among the singular values K r (r = 1 , 2 , . . . , 
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R) is assumed as a matrix Aq obtained by approximating the 

matrix "A" having rank "Q" . 

Q 

A Q =Z K r W r V r 

r=l 

A Q = I» r v I ...(17) 

r=l 

5 Thereby, an approximated vector Xq of the real essential data 
can be estimated using the MP generalized inverse matrix Aq + 
of the approximated matrix Aq by the equation (18) . 
Xq = AQ*b . . . (18) 

As described above, unreliable conditions can be relaxed 
10 by reducing the rank of the matrix "A". However, there are 
some cases where the error is amplified depending on the degree 
of rank reduction. The essential data estimated vector x R 
estimated by the matrix "A" having the rank R which is not 

reduced is expressed by using the equation (13) as follows: 

15 x R = A + b 

= A + (Ax 0 +6) 

X R =S(( X o)r +< ! 6 r )v r ...(19) 

r=l 

As clearly seen from the equation (19) , the smaller the 
singular value, the more the equivalent noise is amplified. 
20 Accordingly, by reducing the rank of the matrix "A" so as 
to discard the smaller singular values, the approximation 
degree for the real value x 0 of the estimated value can be 
increased and the amplification degree of noise can be reduced. 
In the case of using the vector X R obtained by the matrix 
25 "A" of rank "R" as the estimated values of the vector x 0 having 
real values of essential data as components , square estimation 
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error S R is defined as follows by using norm ||x R - Xo|| of -the 
differential vector between the essential data estimated 
vector x R and the real essential data vector x 0 . 

S R =|x R -x 0 f=t((x R ) n -(x 0 ) n ) 2 

n=l 

5 Herein, the equation (19) is substituted. It is assumed that 
1 < r < R < N, and n = 1 to R, R + 1 to N. 

S R =t((x R ) n -(x 0 ) n ) 2 + t((x R ) n -(x 0 ) n ) 2 

n=l n=R+l 

Where, since 1 < r < R in the equation (19) , the vector v R 
= 0 in the case where r > R. Therefore, in the second term 
10 of the above equation which is defined when R + 1 < n < N, 
the vector x R = 0 . Thereby, the equation (20) is obtained. 

S R =£((X 0 ) n + K?B. -(X 0 )J + ±(-(x 0 )J 

n=l n=R+l 

= t«^ 2+ E< X <>)n -(20) 
r=l n=R+l 

The first term of the equation (20) is an error produced 
15 by mixing of noise, and can be estimated by making 8 

substantially constant. On the other hand, it is impossible 
in principle to estimate the second term because the real 
essential data vector x 0 is unknown. However, it is 
qualitatively known that, by reducing the rank W R", the first 
20 term is reduced and the second term is increased. That is, 
there is the optimum rank "R" for estimating the essential 
data vector. 

The optimum value of the rank "R" can be judged as follows . 
The square estimation error S Q between the solution vector 
25 (essential data approximated vector) Xq calculated by the 
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matrix Aq having -the reduced rank and the real essential data 
vector x 0 is expressed by the following equation (21) . 




n=l n=Q+l 

5 Substituting the equation (20) in this equation, 

S Q =S R - ix 2 e*+ ±(x 0 )l ...(21) 

n=Q+l n=Q+l 

From the equation (21) , in the case where S R - S Q is positive, 
the essential data approximated vector xq is judged to have 
a superior estimated value to the original estimated vector 
10 x R , that is, to have a value closer to the real essential data. 

Here, in order to judge whether S R > S Q holds on average 
or not, a probability density function is introduced. 

AQ = (S R )-(S Q ) 
15 AQ= £(K- 2 (e n ) 2 -((x 0 ) n } 2 ) ...<22) 

n=Q+l 

Where < • > represents a probability density function. 
Assuming that the probability density function is represented 
by f (n) = < • >, the probability density function has the 
following character . 

20 £f(n) = l 

n=l 

In the equation (22) , it may be judged that the essential 
data approximated vector xq is a superior estimated value to 
the original estimated vector x R in the case where AQ > 0 . 
In the equation (22) , the probability density functions of 
25 the equivalent noise and the essential data vector are used. 
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In the case of rigorously evaluating AQ, it is desired that 
these probability density functions have been obtained. 

On the other hand, a problem that calculated amounts 
become larger arises for obtaining the probability density 
5 functions for evaluating AQ. Accordingly, in the actual 
calculation, it is desired that an appropriate rank "Q" is 
set by judging whether the solution is good or not while 
adjusting the threshold for discarding the singular values . 
Since this threshold is determined at some degree depending 

10 on the part as an object of ultrasonic imaging and transmission 
and reception conditions of ultrasonic waves , the threshold 
setting unit 41 may be arranged to set the threshold based 
on the scanning part information inputted from the scanning 
part information input unit 40 as shown in Fig. 8. 

15 As described above, by using the MP generalized inverse 

matrix, the linear equations (4) can be solved regardless 
of the condition of the matrix "A" . 

Referring to Fig . 9 again , in the ultrasonic transmitting 
and receiving apparatus, at steps S9 to Sll, the processing 

20 for displaying ultrasonic images is performed based on the 
essential data that has been calculated at step S22 . The 
operation at steps S9 to Sll is the same as that described 
in the first embodiment of the present invention. 

By the way , even when the matrix "A" is square and regular , 

25 the essential data estimation may be performed by using the 
MP generalized inverse matrix. The reason for this is that, 
in the exact method, the following problems are pointed out 



in the lecture on the wave theory ("Wave Theory II" lecture 
handout, "(9) Solving Method of Simultaneous Equations", 
searched on March 14, 2003, via the Internet, <URL: 
http : //www .mobile . ss . titech . ac . jp/~takada/hadou/>) by 
5 Junichi Takada (Department of Electrical and Electronic 
Engineering, Graduate School of Science and Engineering, 
Tokyo Institute of Technology) . 

In the case where "b" includes error Ab in the linear 
equation Ax = b, the variation Ax of the solution "x" is given 
10 by the following equations . 

A(x + Ax) = b + Ab 

AAx = Ab . . . (23) 
Simultaneously, the following relation holds. 

||A|| ||Ax|| > ||Ab|| 
15 HA" 1 !! ||b|| > ||x|| 

Where || • || represents a norm. Therefore, the following 
relationship holds between the relative error of "b" and the 
relative error of "x" . 

||Ax||/||x|| < ||A|| HA" 1 !! ||Ab||/||b|| 
20 In addition, in the case where there is error AA in the 

matrix "A" , the variation of the solution is given by the 
following equation by using the equation (23) . 

(A + AA) (x + Ax) =b 

AA(x + Ax) + AAx = 0 
25 Simultaneously, the following relationship holds. 

IIA" 1 !! ||AA|| ||x + Ax|| > ||Ax|| 
Therefore, the following relationship holds between the 
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relative error of "A" and the relative error of "x" . 
||Ax||/||x + Ax|| < ||A|| HA" 1 )! ||AA||/||A|| 

Here, assuming that a number of conditions regarding 
"A" is represented by cond(A) = ||A|| ||A _1 ||, the relative error 
5 of "x" becomes cond (A) times larger at the maximum when error 
occurs inboth "A" and "b" . By the way , the number of conditions 
can also be expressed by the ratio of the maximum singular 
value (Kmax) and the minimum singular value (iQni n ) when singular 
value decomposition is performed on the matrix "A" . That is , 

10 COnd (A) = Kaax/Kmin. 

The matrix A is defined as an adverse condition when 
the number of condition cond (A) is large . In such case , Takada 
described that, when obtaining the solution vector x, the 
minimum norm solution should be obtained by not using the 

15 exact method, but using the singular value decomposition. 
That is , in the ultrasonic transmitting and receiving method 
according to this embodiment, it is also desired that the 
resolution is obtained by using the MP generalized inverse 
matrix when the matrix "A" is an adverse condition. 

20 Next, an ultrasonic transmitting and receiving method 

according to the third embodiment of the present invention 
will be described by referring to Figs. 10A and 10B. The 
ultrasonic transmitting and receiving method according to 
this embodiment can be used in both of the ultrasonic 

25 transmitting and receiving apparatuses as shown in Figs. 1 
and 8 . 

As shown in Fig. 10A, in the ultrasonic transmitting 

38 



and receiving method according to this embodiment, an 
ultrasonic beam TX is transmitted toward single direction, 
and reception focusing processing is performed with respect 
to the detection signals of the received ultrasonic echoes 
5 so as to form three reception focal points A, B , C. Thereby, 
measurement data y A , y B , yc in. three regions can be acquired. 
In addition, as shown in Fig. 10B, by the transmission beam 
TX and reception beams RX1 , RX2 , RX3 , three different intensity 
prof iles with (a A ,p A ,y A ), (a B ,p B ,y B ), (Oc,Pc,Yc) as components , 
10 respectively, are set. 

Accordingly, the simultaneous equations (26) are 
obtained. Where, x A , x B , x c represent the essential data 
(estimated data) in regions A, B, C. 

Ya = ct A -x A + P A x B + y A -x c 
15 Yb = a B x A + PbX b + y B x c 

y c = cXc X A + Pc-x B + y c -x c ... (26) 

By solving the simultaneous equations (26) , the 
essential data x A , x B , x c can be calculated. The calculating 
method of the essential data is the same as that described 

20 in the first or second embodiment of the present invention. 
By the way, as described by referring to Fig. 6, the number 
of equations constituting the simultaneous equations may be 
increased in order to raise the estimation accuracy of the 
essential data. 

25 According to this embodiment, since only single 

ultrasonic beam is transmitted at a time, the effect of the 
side lobes when performing beam transmission can be reduced, 
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and the side lobe component: Included In the ultrasonic echo 
signals can be reduced by the calculation of the essential 
data. In addition, plural reception focal points are formed 
with respect to the received ultrasonic echoes, the Image 
5 data representing plural parts within the object can be 
acquired at high speed with less number of transmission and 
reception of ultrasonic waves . 

An ultrasonic transmitting and receiving method 
according to the fourth embodiment of the present invention 
10 will be described by referring to Figs. 11A and 11B. The 
ultrasonic transmitting and receiving method according to 
this embodiment can be used in both of the ultrasonic 
transmitting and receiving apparatuses as shown in Figs . 1 
and 8 . 

15 As shown in Fig. 11A, in the ultrasonic transmitting 

and receiving method according to this embodiment , an 
ultrasonic beam TX is transmitted toward a single direction, 
and reception focusing processing is performed with respect 
to the received ultrasonic echo signals by using different 

20 reception delay patterns RX1 and RX2 , and thereby, plural 
pieces of measurement data y^i and Ya2 relating to the region 
A are acquired. In addition, as shown in Fig. 11B, by varying 
the reception delay pattern, intensity profiles PR1 and PR2 
corresponding to the different measurement data Yai and y^ 

25 are set. In order to set the different intensity profiles, 
the aperture condition of the ultrasonic transducer array 
10 when transmitting ultrasonic beams and the transmission 
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condition including the transmission delay pattern may be 
varied other than varying the reception condition including 
reception focusing processing. 

According to the method, the simultaneous equations (27) 
are obtained. Where, "x" represents the essential data 
(estimated data) in the region "A" . 

YaI = CtAi-X 

Ya2 = <Xa2-x . . . (26) 

By solving the simultaneous equations (27) , the 
essential data "x" can be calculated. Here, since the 
simultaneous equations (27) have a larger number of equations 
than a number of unknowns "x", the equations have extra 
conditions. On this account, when solving the simultaneous 
equations (27) , a plausible solution is obtained. The 
calculating method of the essential data is the same as that 
described in the first or second embodiment of the present 
invention . 

According to this embodiment, since the essential data 
relating to the predetermined region is calculated by using 
plural different intensity profiles and plural pieces of 
measurement data, the estimation accuracy of the essential 
data can be raised. 

Next, an ultrasonic transmitting and receiving method 
according to the fifth embodiment of the present invention 
will be described by referring to Figs. 5A and 5B. The 
ultrasonic transmitting and receiving method according to 
this embodiment can be used in both of the ultrasonic 
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transmitting and receiving apparatuses as shown in Figs . 1 
and 8. 

As shown in Fig. 5A, the first multi-beam transmission 
is performed toward TA, TB, TC directions, and reception 
5 focusing processing is performed with respect to the detection 
signals of the received ultrasonic echoes so as to form 
reception focal points on the isochronal surface PHI in TA, 
TB , TC , TD , TE directions . By using thus obtained measurement 
data ym , Ybi, Yci/ Ydi/ ysi and intensity profiles corresponding 

10 to the respective measurement data as shown in Fig. 5B, the 
essential data Xai, x B i, x C i, x D i, x E i in the regions A, B, C, 
D, E are calculated. The calculating method of the essential 
data is the same as that described in the first or second 
embodiment of the present invention. 

15 Then, the second multi-beam transmission is performed 

toward directions including TD, TE, and reception focusing 
processing is performed with respect to the detection signals 
of the received ultrasonic echoes so as to form reception 
focal points on the isochronal surface PHI in TA, TB, TC, 

20 TD, TE directions. By using thus obtained measurement data 
Ya2/ Yb2/ Yc2/ Yd2, Ye2 and intensity prof iles corresponding to 
the respective measurement data, the essential data , x B 2 , 
Xc2/ x D2 , x E2 in the regions A, B, C, D, E are calculated. 
Subsequently, an average or an arithmetic mean is 

25 obtained between the essential data x^, x b1 , Xci, x D1 , x Ei 
calculated by performing the first multi-beam transmission 
and the essential data x^ , x B2 , x c2 , x D2 , x E2 calculated by 



performing the second mult:! -beam transmission, thereby the 
essential data x A , x B , Xc, x D , x E are calculated. 

As described above, by performing the multi-beam 
transmission and reception of ultrasonic waves under 
5 different conditions, and averaging the essential data 

calculated under the respective conditions , the measurement 
error can be reduced, and thereby, the essential data with 
high estimation accuracy can be obtained at high speed. 

By the way, ultrasonic beams may be transmitted and 

10 received over plural times without changing the plural 

directions in which multi-beam transmission is performed, 
so that plural kinds of essential data is calculated with 
respect to each of the predetermined regions A, B, C, D, E, 
and averaged. By such method, the measurement error can also 

15 be reduced. 

Next, a modified example of the ultrasonic transmitting 
and receiving methods according to the first to fifth 
embodiments of the present invention will be described by 
referring to Fig. 1 and Figs. 12 to 14. Fig. 12 is a block 

20 diagram showing the constitution of the ultrasonic 

transmitting and receiving apparatus in which the modified 
example of the ultrasonic transmitting and receiving methods 
according to the first to fifth embodiments are used. This 
ultrasonic transmitting and receiving apparatus is formed 

25 by adding an intensity profile correction unit 36 to the 
ultrasonic transmitting and receiving apparatus as shown in 
Fig. 8. The rest of the constitution is the same as that of 
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the ultrasonic transmitting and receiving apparatus as shown 
in Fig . 8 . 

As described in the first embodiment of the present 
invention, in order to estimate the essential data "x", 
5 measurement data and intensity profiles in the regions defined 
by the transmitting and receiving directions of ultrasonic 
beams and isochronal surfaces are used . Therefore , strictly , 
in order to estimate the essential data relating to all of 
the regions within the object, intensity profiles 

10 corresponding to both all directions TA, TB, ... and all 
isochronal surfaces PHI, PH2 , . . . are required. 

Here, ultrasonic waves are scattered by the relative 
action with the reflector in the transmission paths, and 
attenuated. Accordingly, the intensity profile is thought 

15 to change subtly depending on the time (depth) . Therefore, 
in the modified example , the intensity prof ile on the arbitrary 
time surface is corrected while considering the attenuation 
of ultrasonic waves , and the estimation of the essential data 
is performed by using the corrected profile. 

20 Fig. 13 schematically shows the state of ultrasonic beam 

in the region in which the medium within the object changes. 
As shown in Fig. 13, the acoustic characteristic impedance 
of the medium changes from zi to z 2 on the isochronal surface 
PHI . A part of the ultrasonic beam US propagating the medium 

25 within the object is reflected from the isochronal surface 
PHI , and the other part is transmitted through the isochronal 
surface PHI . 
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The degree of attenuation of an ultrasonic beam can be 
estimated from the intensity of the received ultrasonic echo 
signal. For example, when the ultrasonic echo reflected by 
a certain boundary is large, the amount of ultrasonic beam 
5 transmitted through the boundary is small , and the ultrasonic 
beam is estimated as being largely attenuated afterwards . 
In the case where an ultrasonic wave is transmitted through 
the boundary from the medium having acoustic characteristic 
impedance zi toward the medium having acoustic characteristic 
10 impedance z 2 , the relationship between the reflection 

coefficient R and the transmission coefficient T of ultrasonic 
wave intensity is expressed by the equations (28) and (29) . 

R = 101og|(z 2 - zi)/(z 2 + Zl )| ...(28) 

T = 1 - R . . . (29) 
15 Accordingly, based on the intensity of the ultrasonic echo 
signal of the measurement data on a certain isochronal surface 
and the intensity profile on the isochronal surface, the 
intensity profile on the next isochronal surface can be 
estimated. 

20 Fig. 14 is a flowchart for explanation of the correction 

operation of the intensity profile in the modified example. 
As shown in Fig. 14 , the ultrasonic transmitting and receiving 
apparatus acquires measurement data by transmitting and 
receiving ultrasonic beams at steps SI to S7 , and calculates 

25 essential data based on the measurement data and the intensity 
profile at step S22 . 

At step S30, the intensity profile correction unit 36 



estimates the transmission coefficient "T" by assuming that 
the Intensity of the essential data In the Isochronal surface 
obtained immediately before at step S22 is proportional to 
the reflection coefficient "R" , and corrects the intensity 
5 profiles to be used afterwards based on the transmission 
coefficient "T" in the respective regions on the isochronal 
surface PH2 . The corrected intensity profiles are used for 
the essential data calculation in the isochronal surface PH2 
at step S22 . Alternatively, at step S30, the intensity of 

10 the measurement data may be used instead of the essential 
data. In this case, the calculation of the essential data 
at step S22 and the correction of the intensity profiles at 
step S30 can be performed in parallel. Further, at steps S9 
to Sll, the processing for displaying the ultrasonic image 

15 based on the calculated essential data is performed. The 
operation at steps S9 to Sll is the same as described in the 
first or second embodiment of the present invention. 

As described above, by performing correction on the 
preset intensity profiles as the need arises, precise 

20 essential data can be acquired with respect to each region 
included within the object. 

Such intensity profile correction is effective in the 
case where there is a strong reflector within the object. 
However, actually, in many regions where acoustic 

25 characteristic impedance of the medium changes, the 

reflection coefficient is very small, and most of ultrasonic 
waves can be assumed to be transmitted. Accordingly, the same 



intensity profile may always be used for the same direction . 
Alternatively, intensity profile correction may be performed 
only in the case where there is a strong reflector within 
the object. In this case, the threshold value of the 
5 reflection coefficient required for intensity profile 

correction may be set in the intensity profile correction 
unit 36. 

In the above described embodiments , the case of using 
a two-dimensional transducer array in which plural ultrasonic 

10 transducers are arranged in a matrix form has been described. 
However, in the case where various ultrasonic transducer 
arrays such as a linear array in which plural ultrasonic 
transducers are arranged one-di mensionally , a convex array, 
a concave array, an annular array in which they are arranged 

15 coaxially, ultrasonic image data with high accuracy can also 
be similarly obtained by obtaining intensity profiles 
corresponding to the respective arrays . 

As described above , according to the present invention , 
when the image data relating to the regions included in the 

20 object is obtained, since the essential data is calculated 
based on the measurement data obtained by receiving ultrasonic 
echoes and the profiles as to acoustic pressure intensity 
formed within the object, ultrasonic images with high image 
quality and reduced effect of side lobes can be acquired. 

25 
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